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ABSTRACT 



We present spatially resolved near-IR spectroscopic observations of 15 young 
stars. Using a grism spectrometer behind the Keck Interferometer, we obtained 
an angular resolution of a few milli-arcseconds and a spectral resolution of 230, 
enabling probes of both gas and dust in the inner disks surrounding the target 
stars. We find that the angular size of the near-IR emission typically increases 
with wavelength, indicating hot, presumably gaseous material within the dust 
sublimation radius. Our data also clearly indicate Br7 emission arising from hot 
hydrogen gas, and suggest the presence of water vapor and carbon monoxide gas 
in the inner disks of several objects. This gaseous emission is more compact than 
the dust continuum emission in all cases. We construct simple physical models 
of the inner disk and fit them to our data to constrain the spatial distribution 
and temperature of dust and gas emission components. 

Subject headings: stars:pre-main sequence — stars: circumstellar matter — stars: individual (AS 
205, AS 442, DG Tau, MWC 275, MWC 480, MWC 758, MWC 863, MWC 1080, 
HD 141569, HD 142666, HD 144432, RY Tau, RW Aur, V1295 Aql, VV Ser)— 
techniques:spectroscopic — techniques:interferometric 



Introduction 



Protoplanetary disks provide a reservoir of material f rom which planets may form, and 
the abundance and properties of extra-solar planets (e.g.. lMarcy et al.ll2005l ). as well as the 
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architecture of our own solar system, suggest that planets frequently form in or migrate 
through inner regions of protoplanetary disks. Furthermore, the innermost disk regions 
represent the interface between the inwardly accreting disk and the magnetized ce ntral star. 



and it is here that material accretes inward or is launched in winds or outflows (IShu et al. 



1994J ). Knowledge of the distribution of material in the inner disk is therefore crucial for 
understanding the mass assembly and angular momentum evolution of pre-main-sequence 
stars. 



Modeling of spectral energy distributions (e.g. 



Bertout et al. 




1988 


Lada & Adamsl 


. Naiita et al. 


Il996a: 


Blake & Booeert 



2004 : iNajita et al.l 120071 ) provide important insights into the structure of protoplanetary 
disks in the terrestrial planet forming region. However, since these techniques typically 
rely on spectral information as a proxy for spatial information, they require assumptions 
about underlying geometric, temperature, or velocity structure. Recently, the technique of 
spectro-astrometry, which capitalizes on the fact that emission centroids can be measured 
more precisely th an the available angular r esolution, has enabled less ambiguous constraints 



on disk structure (jPontoppidan et al.ll2008l ). However, substantial gaps in our understanding 



remain due to a lack of high angular resolution observations. 

Near-IR interferometry, which synthesizes a large aperture using two or more smaller, 
separated apertures, can achieve orders of magnitude higher angular resolution than con- 
ventional telescopes, and can spatially resolve disk terrestrial regions. For example, the 
Keck Interferometer, which combines the light from the two 10-m Keck apertures over an 
85-m baseline, achieves a resolution approximately an order of magnitude higher than that 
attained with a single aperture. This means that angular scales of a few milli-arcseconds, 
corresponding to a few tenths of an AU at typical distances to nearby star-forming regions, 
can be spatially resolved with near-IR interferometers. 

Most interferornetric measurements to date have probed only inner disk dust (e.g.. 



Millan-Gabet et al.l 120071 . and references therein) , which typically dominates the near-IR 



emission. To distinguish gas and dust, spectrally dispersed observations are required. Very 



(Eisner et al. 


2007a 


Isella et al. 


20081). ' 



flEisner et al.ll2007al : lEisneii 120071 : iMalbet et all 120071 : iTatulh et al.l 120071 : iKraus et all 12008 



20081 ). The observations showed intriguing evidence that gas and dust are not 



distributed uniformly in inner disk regions. Moreover, while observations of star s less massive 



than a few Mp foun d gaseous emission to be more compact than dust emission (lEisner et al. 



2007al : lEisnerl 120071 ). observations of the Br7 emission line around two more luminous stars 
found the gas to be more extended than the dust, presumably because the Br 7 emission 
traces outflows from these young systems (jMalbet et al. 2007 : TatuUi et al. 2007). A larger 



3 



sample is required to further investigate such potential trends, and to constrain the general 
properties of inner disk gas in young stars. 

Here we present spectrally dispersed near-IR interferometry observations of a sample of 
young stars, including four T Tauri stars and 11 Herbig Ae/Be stars. Our data constrain 
the relative spatial and temperature distributions of dust and gas in sub-AU-sized regions 
of the disks around these stars. 



2. Observations and Data Reduction 



2.1. Sample 

We selected a sample of young stars (Table [I]) known to be surrounded by protoplanetary 
disks, all of whic h have been observed previously at near-IR wavelen g ths wi t h long-baseline 
interferometers (Mi 



2003: Monnier et al. 



an-Gabet et al 



2005 



12001 



Akeson et al. 



2005a 



Eisner et al. 2004, 2005, 2007c: Colavita et al. 



3). All targets except one have been previ 



ously spatially resolved in the near-IR. HD 141569, which has not been resolved, is thought to 
possess a disk with a cleared inner region; however we include the source here to investigate 
whether any spectral features are spatially resolved even though the continuum emission is 
not. 

Our sample includes four T Tauri stars, pre-main-sequence analogs of solar-type stars 
like our own sun, and 11 Herbig Ae/Be stars, 2-10 Mq pre-main-sequence stars. The main 
selection criterion in choosing this sample is source brightness (and because Herbig Ae/Be 
stars are typically brighter than T Tauri stars, our sample has more of the former). Our 
experimental setup imposes limiting magnitudes of ~ 7 at near-IR wavelengths and 
V ^ 12 at optical wavelengths. We also require that sources be at zenith angles of less than 
~ 50°, which excludes from our sample an y sources with 6 5-, — 35°. Our sample includes most 
of the sources in the Herbig-Bell catalog (IHerbig fc Belli Il988l ). as well as several additional 
young stars discovered elsewhere, that meet our selection criteria. 



2.2. Observations 

We obtained Keck Interferometer (KI) observations of our sample on UT 2006 Novem- 
ber 12 and UT 2007 July 3. KI is a fringe-tracking long bas eline near-IR Michelson inter- 



ferometer combining light from the two 10-m Keck apertures (IColavita &: Wizinowichl l2003l : 
Colavita et al.|[2003l ). Each of the 10-m apertures is equipped with an adaptive optics (AO) 
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system that corrects phase errors caused by atmospheric turbulence across each telescope 
pupil, and thereby maintains spatial coherence of the light from the source across each aper- 
ture. The AO systems require sources with V magnitudes brighter than ~ 12. Optical 
beam-trains transport the light from each aperture through an underground tunnel to delay 
lines, beam combination optics, and the detector. 

The two outputs of the beam combiner are sent into a dewar that contains a HAWAII 
array. Interferometric fringes are measured by modulating the relative delay of the two input 
beams and then measurin g the modulate d intensity level of the combined beams during four 



"ABCD" detector reads (IColavitalll999l ). The measured intensities in these reads are also 
used to measure atmosphere-induced fringe motions, and a servo loop removes these motions 
to keep the fringe centered near zero phase. 

KI normally measures science data in a spatially filtered wide-band channel, with a 
second output of the beam combiner, dispersed over four channels using a prism, used for 
group delay tracking. For the observations discussed here, we used a mode where the prism 
is replaced with a grism providing an order of magnitude higher dispersion, and we used 
this spectrally di spersed outp u t for ou r science measurements. The grism, whose properties 



are described in lEisner et al.l (l2007bl ). provides a spectral resolution oi R = 230, with 42 
10-nm-wide channels across the i^'-band. To obtain adequate signal-to-noise for the group 
delay measurement with our grism, which has 10 times more pixels and hence 10 times more 
read noise than the prism that is normally used, we require a star brighter than ~ 7. 

In the following sections we discuss the calibration, modeling, and interpretation of our 
V"^ and flux data. We exclude HD 141569 and VV Ser from most of this discussion. HD 
141569 is unresolved across our observing band, and thus we can say only that the source 
is small compared to the fringe spacing at all observed wavelengths. Conversely, we were 
unable to measure fringes for VV Ser (despite measuring strong K band flux), and for this 
object we can say only that the angular size is large compared to the fringe spacing. We 
defer discussion of these sources to § §4.2f[01 



2.3. Calibration 

We measured squared visibilities (V^) for our targets and calibrator stars in each of the 
42 spectral channels provided by the grism. The calibrator stars are main sequence stars, 
with known parallaxes, whose K magnitudes are within 0.5 mags of the target K magnitudes 
(Tabled]). The system visibility (i.e., the point source response of the interferometer) was 
measured using observations of these calibrators, whose angular sizes were estimated by 
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fitting blackbodies to literature photometry. These size estimates are not crucial since the 
calibrators are unresolved (i.e., their angular sizes are much smaller than the interferometric 
fringe spacing) in almost all cases. HD 163955, a calibrator for MWC 275, is mildly resolved; 
we account for this when computing the system visibility. 

We calculated the system visibility appropriate to each target scan by weighting the 
calibr ator data by the in ternal scatter and the temporal and angular proximity to the target 



data (IBoden et al.l 119981 ). For comparison, we also computed the straight average of the V"^ 
for all calibrators used for a given source, and the system visibility for the calibrator observa- 
tions closest in time. These methods all produce results consistent within the measurement 
uncertainties. We adopt the first method in the analysis that follows. 

Source and calibrator data were corrected for standard detection biases as described by 



Colavital (119991 ) and averaged into 5-s blocks. We accounted for a known bias related to flux 
level by applying an empirically determined correctioE0. Calibrated V'^ were then computed 
by dividing the average measured V'^ over 130-s scans (consisting of 5-s sub-blocks) for targets 
by the average system visibility. Uncertainties are given by the quadrature addition of the 
internal scatter in the target data and the uncertainty in the system visibility. We average 
together all of the calibrated data for a given source to produce a single measurement of 
in each spectral channel. The observations of our targets typically spanned < 1 hour, and 
the averaging therefore has a negligible effect on the uv coverage. 

We investigate the uncertainties in our calibration procedure in several ways. First, we 
examine the calibrated V'^ for each scan for one of our targets (V1295 Aql) where we obtained 
several scans. We also compare calibrated V^^ using one calibrator or another. Differences 
between scans and between the two calibrators provide an estimate of the uncertainties. 
Results of this test are shown in Figure [H Across most of the band, channel-to-channel 
uncertainties are a few percent or less. However, there is a large dispersion in the various 
measurements around 2.05 /xm. 

A second probe of the uncertainties is provided by applying our calibrations to two main- 
sequence check stars, HD 167564 and HD 171149 (Figure [2]). These stars were observed 
as calibrators for VV Ser, and are 30-40° away from the other nearest calibrators in our 
dataset. Thus, we expect the uncertainties derived here to be larger than for our target 
stars, which are within 10° of their calibrators. For these two stars, the standard deviation 
of the calibrated V"^ across the bandpass is 1-2%. However, the data for both stars exhibit 
an apparently systematic feature around 2.05 fim, the same spectral region that exhibits 
larger uncertainties in Figure [H 



http: / /msc. caltech.edu/software/KISupport/dataMemos/fluxbias.pdf 
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We do not have a simple explanation for these large errors at the short- wavelength end of 
our band. Telluric CO2 features lead to absorption (and hence lower photon counts) on either 
side of 2.05 fim, but not at 2.05 /xm; a plot of the atmospheric transmission is shown in Figure 
[31 Furthermore, it is difficult to imagine how any atmospheric or instrumental absorption 
could vary as quickly as the observed variations (and Figured] shows that fluxes in this region 
vary both positively and negatively, contrary to expectations for time- variable absorption). 
We speculate that the observed variability near 2.05 fim may be due to constructive and 
destructive interference caused by a known phase irregularity at this wavelength in a dichroic 
optic. 

We adopt 3% channel-to-channel uncertainties for our target stars. However, Figures 
[T] and [2] indicate that data shortward of ~ 2.05 fim likely has larger errors, and should 
be treated with additional caution. The normalization of V'^ versus wavelength (i.e., the 
average value of V"^ across the band) has an additional uncertainty. Observations of binary 
stars with known orbits show that the calibrated V"^ have a systematic uncertainty of < 3%g. 
We therefore assume that in addition to the 3% channel-to-channel uncertainties described 
above, the normalization of V'^ is uncertain by ~ 3%. 

Calibrated for our sample are shown in Figure |H As discussed above, HD 141569 
and VV Ser are excluded from these plots (and from our subsequent analysis) because they 
are unresolved and over-resolved, respectively, in our observations. 



2.4. Flux Calibration 

We used t he count rates in each channel observed during "foreground integrations'' 



(IColavital Il999l ) to recover crude spectra for our targets. We divided the measured flux 
versus wavelength for our targets by the observed fluxes from the calibrator stars, using 
calibrator scans nearest in time to given target scans, and then multiplied the results by 
template spectra suitable for the spectral types of the calibrators. 

We perform the same tests of our calibration procedure as employed in §2.3[ Variations 
in calibrated fluxes for several scans, using different calibrators, suggest that there are un- 
certainties in the overall slopes of the spectra that lead to channel-to-channel uncertainties 
of ~ 5-10% (Figured]). These slope variations may arise because of different coupling effi- 
ciencies of light at different wavelengths into the single-mode fiber feeding the detector; the 
relative couplings could also change with the atmospheric seeing. Tests of this calibration 



^http: / / msc.caltech.edu/software /KISupport / dataMemos /fl uxbias.pdf 
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procedure for main sequence stars of known spectral type, calibrated using other calibrator 
stars, indicate channel-to-channel uncertainties of a few percent (Figure [2]), and also show 
evidence for these slope uncertainties. 

Because our measured spectra have these potential slope errors, we use broadband pho- 



near-IR flux measurements from the literature (Mendoza 1966 


. 1968: Glass & Penston 19741: 


Hillenl 


orand et al. 


992; 


Kenvon & Hartmann 


1995: Skrutskie et al. 1996: 


Jensen & Mathieu 


1997 




Malfait et a 


. 1998: Eiroa et al. 12002: 


Koresko 2002: ( 


Dutri et al. 


2003: Prato et al. 


2003 




Eisner et al. 


2004 


), including multiple measurements where available, and then fit- 



ted a straight hne to these data. We scaled our measured KI spectra so that the slope and 
normalization matched those of the fitted lines. With these corrections, the spectra contain 
information about narrow spectral features, but do not contain any original information 
about the overall normalization or slope (this information comes from the broadband lit- 
erature photometry). Moreover, since our targets tend to be photometrically variable, this 
procedure introduces some additional uncertainty in the absolute flux level at the epoch of 
our KI observations. 

We estimate uncertainties in our calibrated, slope-corrected spectra of 5-10%. How- 
ever, Figures [1] and [2] both show evidence for larger errors around 2.05 /im. As for the V'^ 
calibration discussed in §2.3[ the fluxes in this spectral region should be treated with caution. 



2.5. Features of Calibrated V"^ and Fluxes 

Figure H] shows the V"^ and fluxes calibrated with the procedure outlined in § §2.3f[^^ 
The calibrated fluxes typically increase with wavelength across the K band, consistent with 
expectations for warm (< 2500 K) circumstellar emission. The V"^ exhibit different behaviors 
with wavelength, ranging from positive to negative slopes. Interpretation of these V"^ trends 
requires an accounting of the differing angular resolution of the observations as a function of 
wavelength, since the resolution gets coarser at longer wavelengths. It is more straightforward 
to explain the V"^ behavior by first converting the measured into angular sizes; we do this 
below in §3.21 

The fluxes and V'^ do not always behave monotonically across the observing window, 
due to the presence of spectral features associated with hot hydrogen gas, and warm CO and 
H2O vapor. For example, MWC 480, MWC 275, and V1295 Aql show clear bumps in both 
flux and V"^ near 2.165 /im, associated with the Br7 transition of hydrogen. Spectral features 
like these may arise in absorption in stellar photospheres, or in emission in the circumstellar 
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material. In fj3|, we determine contributions to the fluxes and V'^ from the central star, 
enabling investigation of the spectral features arising from the circumstellar emission. 



3. Modeling 

In this section we use our V'^ and flux measurements to constrain the distribution of dust 
and gas within 1 AU of our sample stars. We first consider the simplest possible model for the 
emission, a uniform disk ( §3.21) : this model provides an estimate of how the angular size of the 
emission depends on wavelength, with essentially no model assumptions. We then explore 
more physically realistic models that include both gas and dust components. Specifically, 
we model our sources as optically thin gaseous accretion disks that extend as far as the 
dust sublimation radius, at which point emission from the puffed-up dust sublimation front 
dominates. We include various gaseous species in the inner disk as needed to fit the data. 
Because our V"^ and flux measurements contain contributions from both the circumstellar 
disks and the unresolved central stars of our targets, we remove the stellar component from 
our measurements before fitting the models described above. The procedure by which stellar 
and circumstellar components are separated is described in §3.11 



3.1. Separating Stellar and Circumstellar Components 



Before fitting physical models to the data, we remove the contribution of the central 
stars from our measurements. This enables us to model only the circumstellar material 
around our sources. Because the central stars are unresolved, we know that = 1, and we 
need only the ratio of the stellar and circumstellar fluxes to remove the stellar contribution 
to the measured V'^ and fluxes. 

We estimate the circumstell ar-to-stellar flux ratio at each observed wavelength using 
spectral decomposition (see, e.g.. iMillan-Gabet et al.ll200ll ). Using optical photometry from 
the literature, we fit the stellar photosphere. For each source we assume the spectral type 
given in the literature (also listed in Table [T]) and determi ne the stellar radiu s and r eddening 
providing the best fit to the data; the reddening law of ISteenman fc Thg (Il99ll ) is used. 
In general we assume that the optical {VRI) photometry traces the un-veiled stellar pho- 
tosphere. H o wever , for AS 205 A, we have veiling measurements at R and / bands from 
Eisner et al.l (120051 ) . and we use these to precisely fit the photospheric flux. We then extrap- 
olate the photosphere to the K band using Kurucz models, which include stellar spectral 
features like Br7 absorption in A stars or CO overtone absorption in cooler stars. We com- 
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pare the extrapolated stellar fluxes with our measured i^'-band fluxes to determine stellar 
and circumstellar contributions. 

With these circumstellar-to-stellar flux ratios, we remove the stellar components of the 
visibilities: 



disk - s ^ IT' r ' y^) 

where -Fdisk/-^* is the circumstellar-to-stellar flux ratio. We flt our models for the circum- 
stellar emission to Fdisk and V^^^^. below. 



3.2. Uniform Disk Sizes 



Before interpreting our data in the c ontext of physical models, we begin with a simple 
geometrical model: a uniform disk (e.g., Eisner et al.ll2003l ). We flt the Vj^^ data for each 
source, in each channel, with the model. The results, shown in Figure [SI give the "spectral 
size distribution" of the i^-band emission, illustrating how the spatial scale of the near-IR 
emission depends on wavelength. 

Figure [5] shows that the angular diameter of the near-IR emission appears to increase 
with wavelength for all sources in the sample (although the trend is marginal for DG Tau 



and RW Aur, where the data have low er signal-to-noise). As in previous work (lEisner et al. 



2007al : lEisnerl 120071 : iKraus et al.ll2008l ). we interpret this slope as evidence for a hot emis- 
sion component interior to the inner edge of the dust disk. This hot, compact component 
increasingly dominates the observed emission at shorter wavelengths, and hence the angular 
diameter of the emission appears to decrease. 

In addition to showing increasing size with wavelength. Figure [5] shows that several 
objects exhibit emission from the Br7 transition of hydrogen that is more compact than the 
circumstellar continuum em ission. We s ee indications of compact Br7 emission from MWC 
480 (reported previously by lEisnerll2007l ). MWC 758 (marginally), HD 144432 (marginally), 
MWC 863, MWC 275, V1295 Aql, and MWC 1080. 

Finally, we see evidence of angular diameter changes due to CO opacity in RW Aur. Fig- 
ure [6] shows a zoomed- in view of the spectral region encompassing several CO ro- vibrational 
overtone {Av = 2) bandheads. The flux increases in each of the CO bandheads, demonstrat- 
ing the presence of CO emission in this object. Moreover, the angular diameter appears to 
decrease within each of the bandheads, indicating that the CO emission is more compact 
than the continuum. 
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3.3. Accretion Disk Models 



As a more physically realistic model, we begin with a dusty circumstellar disk whose 
emission is dominated by the hottest dust near the sublimation radius, and then include 
gaseous emission interior to the inner edge of the dust. We consider several different gaseous 
opacity sources, including continuum emission, Br7 emission, water vapor, and CO. For DG 
Tau and RW Aur, we do not consider V"^ data measured shortward of 2.1 fim, since low 
fluxes (< 1 Jy) lead to extremely noisy data at these wavelengths (Figure H]). 

We do not consider models that include only gaseous emission because previous in- 
vestigators have shown that ring-li ke dust sublimation fronts are needed to fit broadband 
SEDs and in t erferom etry data (e.g.. iMuzeroUe et al.ll2003l : lEisner et al.ll2004j ). Furthermore, 



Eisner et al.l (l2007al ) showed that gas-only disk models, while they can potentially reproduce 
observed spectrally dispersed visibilities, can not simultaneously reproduce observed SEDs. 
We therefore do not consider gas-only models here. 

We do not fit the broadband SED simultaneously with our 2.0-2.4 /im dataset here, to 
avoid additional complication. Since our modeling is not constrained by SEDs outside of the 
K band, the inferred properties of the circumstellar emission may not perfectly refiect reality. 
This is especially true for the gaseous emission, which is often more compact than the spatial 
resolution of our observations. Nevertheless, the modeling presented below demonstrates the 
presence of gaseous emission components, and provides a rough estimate of the size scales 
and temperatures of these components. 



3.3.1. Dust Continuum Emission 



We assume that the dust emission in our model arises from a single-temperature ring 
of matter at the disk radius where temperature s become hot enough for dust sublimation 



(as in previous stu dies; e.g., lEisner et al.ll2004l ). This ring of emission approximates well 



(lEisner et al.l |200J) the pu ffed-up inner edge expe cted for directly irradiated inner disks 
in hydrostatic equilibrium (jDuUemond et al" 2001 ). The width of the ring is assumed to 
be two tenths of its r adius, consistent wit h calculations of the expected width of the dust 
sublimation front (e.g., Isella fc Natta|[2005 ). The free parameters of this model are the radius 
and temperature of the emission ring, -Rdust and Tdust- We determine the best-fit values of 
these parameters using a grid-based y'^ mini mization. Uncertainties are determined directly 
from the surface (e.g.. lEisner et al.ll2004l ). 
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3.3.2. Gaseous Continuum Emission 



We assume the radial temperature and surface density profiles of the gas in the inner 
disk are described by power laws. We take T{ R) oc as appropriate for optically-thin 



disk material (e.g.. IChiang fc GoldreichI 119971 ): T[R) is the radial temperature profile and 



R is the stellocentric radius. Details of the gaseous emission mechanism may alter this 
relationship, but we will consider the radial temperature law fixed for simplicity. We assume 
that the surface densit y is described by S(_R) oc i?"^/^, as inferre d for the protosolar nebula 
(I Weidenschilling] 1 19 771') and assumed in other disk models (e.g., IChiang fc GoldreichI 119971 : 
DuUemond et al.ll200ll ). The optical depth, r, is directly related to the surface density, and 
hence t^{R) oc Depending on the normalization of the optical depth profile, the gas 

may be entirely optically thin, entirely optically thick, or thick in the innermost regions only. 

We define the radial temperature and optical depth profiles over an inner disk radius, 
i?in, to an outer disk radius, -Rout- We take i?out to be the radius of the dust ring described 
in §3.3.11 i?in is left as a free parameter. 

The free parameters of the gaseous disk model are the inner radius, i?in, the temperature 
at the inner radius. Tin, and the optical depth at the inner radius. Tin. When combined 
with the dust component, the dust+gas continuum emission model has five free parameters. 
We determine the values of these parameters that provide the best-fit to the data using a 
Marquardt-Levenberg algorithm. Uncertainties on best-fit parameters are computed from 
the covariance matrix. 

If the gas is optically thick, then only the radial temperature profile is important in 
determining the resultant flux. However for most sources, simple disk models with optically 
thick gas predict fluxes larger by an order of magnitude than those observed from our sample. 
We therefore expect the gas to be (at least partially) optically thin in our modeling. The 2 fim 
emission from optically thin gaseous disk models arises predominantly from the innermost 
radii. Models incl uding dust and opt ically thin gas thus resemble two-ring models (like 
those considered in lEisner et al.l l2007al ) , where hot gaseous emission appears in a ring near 
the inner radius and warm dust emission arises from the ring-like dust sublimation front. 



3.3.3. Water Vapor Emission 

Several objects in our sample show additional flux, and accompanying broad features 
in V"^, at wavelengths < 2.1 /xm. These features are (in some cases) broader than the 
systematic features discussed in §2.21 although the data is still affected by these short- 
wavelength systematic effects. In a previous paper, we reported such features for MWC 480, 
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and modeled the data by addin g a single-tem perature ring with the opacity of water vapor 
to a two-ring dust-|-gas model (lEisnerl 120071 ). Here we refine this simple model, and test 
whether the data for our sources can be fitted with a model where water vapor emission 
occurs over a range of radii from an inner radius where water is thermally dissociated out to 
the dust sublimation radius. 



We use opacities measured for water vapor (lLudwig||l97ll ) at a range of temperatures (~ 
1000-3000 K) appropriate to protoplanetary disks. Using the gaseous disk model described 
in §3.3.2[ we determine the optical depth at each radius using the assumed surface density 
profile, which scales as oc t,AR ) = T,(R)k,u. For disk radii where T < 3000 



K, water vapor can exist and we use k^, from iLudwid (Il97ll ). At smaller stellocentric radii, 
where water would be dissociated, we retain the assumption of continuum emission. While 
it is also possible that water emission, where is exists, adds to rather than replaces the 
underlying continuum, we do not consider this possibility in our simple model. 



3.3.4- Brackett Gamma Emission 

Figures H] and [5] show that Br7 emission signatures are present in the fiux and data 
for several objects. In fact, fiat spectra in the vicinity of the Br7 wavelength (2.165 /im) 
may indicate the presence of Br7 emission for some objects. Since more massive stars, 
including essentially all of the Herbig Ae/Be stars in our sample, have photospheric Br7 
absorption, non-detection of this absorption feature in the spectra implies that circumstellar 
Br7 emission (and continuum excess) is filling in the feature. The effects of photospheric Br7 
absorption are included in the Kurucz models we used to separate stellar and circumstellar 
components of the measured V"^ and fiuxes ( §3.11) . 

We model circumstellar Br7 emission by including additional fiux in the gaseous com- 
ponent of our models at 2.165 fim. Since Br7 traces very hot gas (~ 10^ K), we assume 
that the emission arises from the innermost disk annulus. In fact, Br7 emission may be even 
more compact, but we can not constrain such small size scales with the available angular 
resolution. Inclusion of Br7 emission in the model introduces as an additional free parameter 
the gaseous line-to-continuum ratio in the innermost gaseous annulus. 
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Results 



4.1. Dust and Gas within 1 AU of young stars 

Our modeling shows clearly that single-temperature rings of emission can not fit the data 
for our sample (Figure [71 Ta ble El) . This confirms preyious inferences from spectrally dis- 



perse d interferometric results (jEisner et al.ll2007al : lEisnerll2007l : iKraus et al 



2008 



Isella et al. 



20081 ) and from high angular resolution continuum observations (ITannirkulam et al.ll2008l ). 
The implication is that models where all of the near-IR emission arises from the inner edge 
of the dusty disk are untenable. Rather, a warm dusty component and a hotter, presumably 
gaseous, component at smaller stellocentric radii are required to fit the data. 

The simple dust+gas continuum model ( §3.3p predicts V'^ and fluxes compatible with 
the observations for most sources (Figure [71). We also consider a model where the gaseous 
opacity is due to water vapor at stellocentric radii where H2O can potentially exist. As 
shown in Figure U\ these models generally produce fits of comparable quality to the simple 
continuum models. Models including water vapor also produce more physically realistic 
fitted values of -Rdust:- the fitted dust temperatures are between ~ 1200 and 1500 K, in 
agre ement with expected sublimation temperatures for silicate dust in protoplanetary disks 
(e.g. [Pollack et al.lll994l ). However, for RY Tau, the simple dust+gas (continuum) model 



produces a superior fit to the data. 

We illustrate the relative contributions of gas and dust to our models in Figure [HI The 
figure shows the fluxes produced by the gaseous component and the dust component for our 
best-fit models including water vapor opacity. The gaseous component includes continuum 
emission, water vapor emission, and Br7 emission. In all cases, the gaseous flux represents 
a significant fraction of the total. 

For all sources except AS 205 A, the best-fit dust+water model requires emission from 
material substantially hotter than 3000 K (Table [21), implying the presence of substantial 
gaseous continuum emission. The inferred inner disk temperatures for models including 
water vapor opacity are, in most cases, higher than for models assuming only gas continuum 
opacity. This difference suggests that the best-fit for models including water vapor requires 
the gaseous emission to be dominated by the hot continuum component rather than the 
cooler water vapor component. For AS 205 A, in contrast, it appears that water vapor alone 
can explain the observed gaseous emission from the inner disk. 

RW Aur shows evidence of spatially resolved CO in the inner disk. Emission from the 
At; = 2-0, 3-1, and 4-2 bandheads is detected in our spectra, and appears to be more 
compact than the surrounding continuum emission (Figure [HJ. This suggests an origin of 
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the CO interior to the dust sublimation radius. The inferred size of the CO emission for 
this object is compatible with th at inferred from mo deling of a high dispersion spectrum in 
terms of a Keplerian disk model (iNajita et al.ll2003l ). 



DG Tau, which is also known to possess CO overtone emission (e.g., ICariill989l ). shows 
no evidence of such emission in our data. This discrepancy may be due to the poor quality 
of our data for this object or to variability. The CO eraission from DG Tau is k nown to be 
variable from studies at multiple epochs fe.g.. lCarrlll989l:lNaiita et al.ll2000U2003l ). and some 
previous studies failed to detect the emission (e.g.. ICreene &: Lada 1996 ). 



MWC 758, MWC 275, V1295 Aql, and AS 442 all show excess flux long-ward of ~ 2.3 
/im, consistent with the presence of CO overtone emission (Figure H]). However, none of these 
show evidence in the visibilities of size differences between CO emission and continuum. 
Given the signal-to-noise and spectral resolution of our current data, we do not attempt to 
model the CO emission from any of these targets here. 

Several objects exhibit Br7 emission from hot hydrogen gas, and this emission is more 
compact than the surrounding continuum emission in all cases. Table [3] lists the inferred 
stellocentric radii from which the Br7 emission originates. However, we have somewhat 
arbitrarily placed the Br7 emission at the innermost edge of the gaseous accretion disk in 
our models ( §3.3.41) . Since the angular resolution of our observations is ~ 0.1 AU, and there 
is some degeneracy in our models between the Br7 flux and angular scale, we can only state 
confidently that the Br7 emission arises from radii less than 0.1 AU. 

While our model fits indicate that Br7 emission is present around RW Aur, and is more 
compact than the continuum, we are not confident in this result given the large uncertainties 
on V"^ for this source in this spectral region. The data for DG Tau also show evidence 
for compact Br7 emission, but the noisy data in this spectral region argue for a cautious 
interpretation. Both DG Tau and RW Aur have previousl y detected (spatially unresolved) 
Br7 emission ( iFolha fc Emerson! l200ll : iNajita et al.l Il996bl ) , providing some support for our 
tentative detections. 

MWC 480, HD 144432, MWC 863, V1295 Aql, and MWC 275, all targets that ex- 
hibit compact Br7 emission, were previously reported to sho w strong Br7 emission in 



high-dispersion (spatially unresolve d) spectroscopic observations (iGarcia Lopez et al.l 12006 



Brittain et al.ll2007l : lBerthoudll2008h. RY Tau. MWC 758, HD 142666, and AS 205 also have 
previously reported B r7 emission ( Najita et al.lll99^ : Folha fc Emerson 2001 : Garcia Lopez et al 
20061 : lBerthoudll2008l ) . However, the equivalent widths of the Br7 emission are small relative 
to other sources (e.g., MWC 275 or DG Tau), and it is not surprising that we do not detect 
(or only marginally detect) Br7 emission from these objects here. 
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4.2. An unresolved source: HD 141569 



Previous near-IR interferometric observations of HD 141569 found it to be unresolved 
in the K band, suggesting that the K band emission is produced entirely by the unresolved 
central star. We re-observed it here to search for spatially extended circumstellar emission 
associated with gaseous spectral features, which might have been washed out in the previous, 
broadband observations. The observations presented here, however, show that the source 
is unresolved at all wavelengths. We also see no evidence of flux above the level expected 
from the stellar photosphere at any observed wavelength. We conclude that the emission 
in each of the spectral channels included in our observations is more compact than ~ 0.1 
AU (at the target distance), as expected for stellar emission. We may also infer that no 
spectral channels contain emission brighter than ~ 1% of the stellar flux within the 50 mas 
field of view of KI (corresponding to stellocentric radii of ~ 2.5 AU for this target), since the 
incoherent contribution of any stronger extended emission would have reduced the measured 
visibilities. Our finding is cons istent with previous studies that foun d gaseous emission only 
at stellocentric radii > 10 AU (iGoto et al.ll2006l : iBrittain et al.ll2007l ). with no warm gas that 
could have been detected in our observations. 



4.3. An over-resolved source: VV Ser 



One of the objects in our sample, VV Ser, appears to be over-resolved in our observa- 
tions. That is, while we detect flux from the object, its angular size is large compared to 
the ~ 5 mas fringe spacing and hence its is un-detectably close to zero. Because we can 
only place a lower limit on the size of the emission, and can not trace how the size depends 
on wavelength, we have excluded VV Ser from the analysis presented above. 

Previous observations of VV Ser, with fringe spacings comparable to those obtained in 
the present study, measu red a uniform disk angu lar diameter of ~ 4 mas, and showed the 
disk to be nearly edge-on (lEisner et al.ll2003l . 120041 ). Our present non-detection indicates that 
the size scale of the K band emission has varied signiflcantly between 2003 and 2007. The 
sourc e is also known the photometrically variable at optical through infrared wavelengths 
(e.g., 



Herbst fc Shevchenkolll999l : Eiroa et al.ll2002[ ) 



These flndings may indicate variability associated with the UX Ori phenomenon. In this 
scenario, a vertically extended inn er edge of the nearly edge-on disk around VV Ser period- 
ically blocks the central star (e.g., iDuUemond et al.ll2003l : iPontoppidan et al.l 120071 ). When 
the star is blocked, the near-IR emission would be dominated by the extended circumstellar 
component, whereas when the star is visible, the observed size would be the flux- weighted 
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average of the unresolved star and the disk. 



Discussion 



5.1. Compact Br7 Emission 



We observe Br7 emission from several of our targets, and in all cases this line emission 
appears more compact than the surrounding c ontinuum. This con trasts with previous results 
that foui id Br7 emission mo re extended than (IMalbet et al.ll2007l ) or on a comparable spatial 
scale to (ITatuUi et al.ll2007l) the dust co ntinuum. Sources where Br7 appears more extended 
are high-mass stars, and lEisnerl (120071 ) speculated that the Br7 emission from young stars 
may trace both infalling and outflowing components, with the latter increasingly dominating 
for higher mass stars. 

The results presented here belie this hypothesis. We observe Br7 emission more compact 
than dust continuum in sources spanning a range of spectral types, from K3 to BO. Based 
on our current findings, we suggest that young stars typically produce most of their Br7 
emission close to their central stars in accretion columns and/or shocks. Br7 emission in 
extended winds seems rarely to be strong enough to dominate the line emission. However, 
given the small number of Herbig Be stars observed to date, the picture is less clear for 
higher-mass (> 10 Mq) stars. 



5.2. Inner Disk Gas: Trends with Stellar Luminosity 

For the majority of our sources, models that include only dust emission do not fit the 
data well, while models including hot gaseous emission interior to the dust sublimation do 
fit the data ( §4.ip . However, if we examine only the less luminous stars in our sample, the 
T Tauri stars, we find that this conclusion is less robust. For all of the T Tauri stars in our 
sample (RY Tau, DG Tau, RW Aur, and AS 205), dust-only models can provide relatively 
good matches to the slopes seen for fluxes and V^. Models including gas and dust are still 
superior, but the difference is not nearly as pronounced as for the more massive Herbig 
Ae/Be stars in our sample. 

We suggest a simple explanation for this trend. For T Tauri sources, the less luminous 
central stars provide less heating and hence the dust can exist closer to the star than for 
the more massive Herbig stars. While gaseous disks are likely to extend in to similar radii 
in both types of sources, the gas may be hotter near to more luminous stars. Higher mass 
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stars have a larger temperature difference between hot gas and warm dust, and also a larger 
spatial separation of the two components, producing larger temperature gradients. These 
larger gradients, in turn, lead to larger slopes in the and fluxes versus wavelength. 

This is not to say that T Tauri stars do not have significant gaseous emission. Models in- 
cluding gaseous emission provide superior fits to our data, and previous author s have argued 



that g aseous emission is needed to self-consistently model SEDs and visibilities (lAkeson et al. 



2005bl ). Furthermore, modeling of observed gaseous emission line profiles from T Tauri 
stars under the assumption of Keplerian rotatio n has provided ev idence for gaseous emis- 



sion interior to the dust sublimation radius (e.g.. iNajita et al.ll2007l ). However, the observed 
temperature gradients between dust and gas appear less pronounced than do those around 
higher-mass stars. 



5.3. The Nature of the Gaseous Opacity 



Water is predicted to be abundant in the inner regions of disks (e.g.. lGorti fc 



ioUenbach 



2004 ). although it can be red uced in abundance in highly irradiated disks (e.g., iThi &: Bik 
20051 ) . iMuzeroUe et al.l (120041 ) predict that dust-free inner disks will coi nmonly show wate r 
in emission or absorption depending on the disk accretion rate (see also lCalvet et al.lll99ll ). 
However, near-IR water spectral fe atures are detected only rarely in T Tauri stars (e.g., 
Carr et al.l l2004l: Naiita et al.l 120071 ) and have only been reported in a few Herbig Ae/Be 
stars JThi fc Biklboosl . Naiita et al. 2008). 



For most sources, a model that includes water vapor can fit the data well. For AS 205 
A, it appears that water vapor alone can reproduce the observed gaseous emission interior to 
the dust sublimation radius. This is consistent with the detection of abundant water vapor, 
and an inferred ori gin at stellocentr ic radii as small as 0.3 AU, in high dispersion L-band 
spectra of AS 205 JSalvk et ahlboosh . 



However, in some cases (e.g., for RY Tau), models including water vapor provide fits of 
lesser quality than models where the gas emits only continuum radiation. Furthermore, for 
all sources except AS 205 A, our best-fit models require continuum emission in addition to 
water vapor emission, since the gas is inferred to be hotter than 3000 K in the innermost 
regions. This casts some doubt as to whether H2O emission is needed to fit the data in all 
cases. 

For the columns of water vapor implied by our best-fit models, one might expect our 
targets to show H2O emission lines in lower spatial resolution spectroscopic data. High 
dispersion spectroscopic observations of several of our sample objects in the i^-band failed 
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to detect strong water emission lines (Najita et al. 2008; iMandell et al.ll2008l ). and archival 
ISO/SWS data for many of our targets do not show longer- wavelength H2O features that 
typically accompany the 2 fim features modeled in this work. These null results provide 
further reason to question whether water vapor is a viable explanation for the observed 
compact, hot, circumstellar continuum emission from some of our sample. 

We therefore pose the question: what can produce gaseous continuum emission in the 
inner regions of protoplanetary disks? We discuss several possible explanations below, and 
argue that emission from free-free transitions of Hydrogen and from the negative Hydrogen 
ion is the most viable mechanism. 



5.3.1. Refractory Dust Grains 



While dust grains produce continuum emission, the inferred temperatures of the hot 
component (> 3000 K) in our data ar e higher than the s ublimation temperatures for even 



the most refractory dust grains (e.g.. iPoUack et al.lll994l ). Even for calcium, magnesium. 



or titanium rich oxides, sublimation t emperatures above 2000 K occur only for ambient 



pressures higher than 0.1-1 bar (e.g., iLewid 119971 ). orders of magnitude larger than the 



pressures expected in protop lanetary disks (< 0.01 bar for typical densities and temperatures 



in inner disk regions; e.g., iMuzeroUe et al.l 120041 ). Dust emission therefore seems highly 



implausible as an explanation for the hot, compact emission. 



5.3.2. High—n Atomic Hydrogen 

Recombination of photoionized hydrogen into high-n states may also produce a (pseudo) 
continuum opacity. This mechanism is also probably responsible for the Br7 emission ob- 
served toward many of our targets, since it produces hydrogen in the n = 7 level, which 
can then cascade down to the n = 4 level and produce the line emission. Even higher— n 
states, which are longer-lived, might be able to produce more extended and continuum-like 
emission. Photoionization cross sections for these high-n states are ~ 10~^^ cm~^, and so 
assuming a gas column density of ~ lO^'' cm~^ (a typica l value for young stars accreting 



material at ~ 10 ^ Mq yr ^; e.g., iMuzeroUe et al.ll2004l ). we see that a fractional abun- 



dance of these high—n states of > 10~^ would be needed to produce an optical depth larger 
than unity. Estimating the fractional abundance requires knowledge of the ionization and 
radiation structure of the disk, which is beyond the scope of this paper. However, it seems 
plausible that recombination to high—n hydrogen states could contribute some opacity in 



- 19 - 



the inner disk. 



5. 3. 3. Free-free 

Free-free radiation, produced in ionized winds or accretion columns, or in disks (as for 
classical Be stars), may also produce continuum emission. Free-free emission may occur 
around the higher-mass stars in our sample, which produce more intense ionizing radiation 
fields, or in shocked, infalling gas around lower-mass stars. The cross section of free-free 
emission at wavelengths around 2 /im is ~ 10~^^ cm^ dyne~^ at temperatures of a few 
thousand Kelvin. Assuming the free-free emission arises in a disk with a gas temperatuture 
of 3000 K, a gas column density of 10^^ cm~^, a fractional ionization of 10~^, and H/H2=0.01, 
we find an electron pressure of ~ 0.01 dyne cm~^, and an optical depth of ~ 0.01. Fractional 
ionizations of 10~^ can arise from the ionization of metals with low ionization potentials, 
such as Na or Fe; higher fractional ionizations are difficult to achieve since photons with 
> 11 eV are unlikely to penetrate deeply into such a dense disk. However, somewhat higher 
optical depths of free-free emission may be achieved if the gas column density is higher 
than the value assumed above. Some of our targets may be accreting faster than 10~^ 
yr~^, leading to higher gas columns and more free-free emission. It thus seems possible that 
free-free emission could contribute significantly in the inner disk. 



5.3.4. H- 



The negative hydrogen ion, H , appears as a promising candidate for explaining this 



emiss i on, since it can produce f r ee-free cont i nuum emiss ion in the K band (e.g.. lChandrasekhar &: Breen 
19461 : iBell fc Berringtoru 119871 : iJohnI Il988l : iGrayl Il992l ). For free-free emission from H~ at 



3000 K in the i^'-band, the optical depth is (following lGraylll992l ) 



r 



10 



-38 



Up 



cm 



N, 



H 



cm 



-2 



(2) 



With the same assumptions at in §5.3.3[ we find an optical depth of ~ 0.01 in the inner disk. 
As above, higher gas columns may lead to correspondingly higher optical depths. Free-free 
emission from H~ therefore seems like a plausible explanation for the hot, compact emission 
seen in our data. 
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6. Conclusions 

We presented spatially resolved near-IR spectroscopic observations that probed the gas 
and dust in the inner disks around 15 young stars. One source, HD 141569, was unresolved 
at all wavelengths between 2.0 and 2.4 /xni, indicating a lack of dust or gas in the inner 
disk regions. Another target, VV Ser, was over-resolved, indicating that the emission spans 
angles larger than ~ 5 mas at all observed wavelengths. 

The near-IR emission from the remaining targets was resolved, and our data show that 
the angular size of the near-IR emission increases with wavelength in all cases. This behavior 
suggests temperature gradients in these inner disks, arising from the combination of warm 
dust at its sublimation temperature and hotter, presumably gaseous material within the dust 
sublimation radius. Our data clearly indicate emission from the Br7 transition of hydrogen 
in several objects, and suggest that water vapor and carbon monoxide gas are present in the 
inner disks of some targets. 

We constructed simple physical models of the inner disk, including dust and gas emis- 
sion, and we fitted them to our data to constrain the spatial distribution and temperature of 
dust and gas emission components. We considered models including only dust emission; dust, 
gas continuum, and Br7 emission; and dust, gas continuum, water vapor, and Br7 emission. 
Models incorporating only dust emission can not fit the data for any of our sources well. In 
contrast, models including dust and gas emission are suitable for explaining our data. The 
inclusion or exclusion of water vapor in these dust-|-gas models did not substantially affect 
the quality of the fits in most cases. 

For all sources where Br7 emission is observed, we find it to be compact relative to the 
continuum emission. This contrasts with previous findings, which found Br7 emission to be 
extended relative to the continuum around some high-mass stars. The results presented here 
suggest that Br7 commonly traces infalling material around young stars spanning a large 
range in stellar mass. 

CO emission is tentatively observed towards several objects, and we see evidence that 
this emission has a more compact spatial distribution than the dust around RW Aur. For 
other objects, our data are insufficient to place meaningful constraints on the relative spatial 
distribution of CO and other emission components. We will re-observe these targets in the 
near future with higher dispersion, to obtain better signal-to-noise for the relatively narrow 
CO lines and better constrain their spatial distribution. 

While models including water vapor opacity often fit our data well, the best-fit models 
generally also require continuum emission from material that is too hot to be water (since 
water dissociates at ~ 3000 K). We do not have a ready explanation for the source of this 
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hot continuum emission, but we speculate that it may trace free-free emission from hydro den 
and/or H~. The gas densities and fractional ionizations required to produce such emission 
seem plausible in the inner regions of protoplanetary disks, suggesting that free-free emission 
from H and H~ is a viable explanation for the compact continuum emission seen in our data. 

Data presented herein were obtained at the W. M. Keck Observatory, in part from 
telescope time allocated to the National Aeronautics and Space Administration through the 
agency's scientific partnership with the California Institute of Technology and the University 
of California. The Observatory was made possible by the generous financial support of the 
W. M. Keck Foundation. The authors wish to recognize and acknowledge the cultural role 
and reverence that the summit of Mauna Kea has always had within the indigenous Hawaiian 
community. We are most fortunate to have the opportunity to conduct observations from this 
mountain. This work has used software from the Michelson Science Center at the California 
Institute of Technology. The authors thank the entire Keck Interferometer team for making 
these observations possible. We also wish to thank the referee, Geoff Blake, for his thoughtful 
and detailed referee report, which greatly improved the manuscript. 
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Table 1. Target and Calibrator Properties 



Source 


a 


5 


d 


Spectral Type 


mv 


ruK 


References 




Target Stars 


RY Tau 


04 21 57.409 


+28 26 35.56 


140 


Kl 


10.2 


5.4 


1 




DG Tau 


04 27 04.700 


+26 06 16.20 


140 


K3 


12.4 


7.0 


2 




MWC 480 


04 58 46.266 


+29 50 37.00 


140 


A2 


7.7 


5.5 


3 




RW Aur A 


05 07 49.568 


+30 24 05.161 


140 


K2 


10.5 


7.0 


2 




MWC 758 


05 30 27.530 


+25 19 57.08 


140 


A3 


8.3 


5.8 


3 




HD 141569 


15 49 57.75 


-03 55 16.4 


99 


B9/A0 


7.1 


6.7 


3 




HD 142666 


15 56 40.023 


-22 01 40.01 


116 


A8 


8.8 


6.1 


4 




rlU 1444-32 


16 06 57.957 


-27 43 09.81 


145 


A n 

A9 


8.2 


5.9 


4 




AS 205 A 


16 11 31.402 


-18 38 24.54 


160 


K5 


12.1 


6.0 


5 




MWC 863 A 


16 40 17.922 


-23 53 45.18 


150 


A2 


8.9 


5.5 


4 




MWC 275 


17 56 21.288 


-21 57 21.88 


122 


Al 


6.9 


4.8 


4 




VV Ser 


18 28 47.860 


+00 08 40.00 


310 


AO 


11.9 


6.3 


3 




V1295 Aql 


20 03 02.510 


+05 44 16.68 


290 


B9/A0 


7.8 


5.9 


3 




AS 442 


20 47 37.470 


-1-43 47 24.90 


600 


B8 


10.9 


6.6 


3 




MWC 1080 


23 17 25.574 


+60 50 43.34 


1000 


BO 


11.6 


4.7 


3 








Calibrator Stars 








Applied to: 




HD 27777 


04 24 29.155 


+34 07 50.73 


187 


B8V 


5.7 


6.0 


RY Tau.MWC 480,MWC 758 


HD23642 


03 47 29.453 


+24 17 18.04 


110 


AOV 


6.8 


6.8 


DG Tau,RW Aur 




HD23632 


03 47 20.969 


+23 48 12.05 


120 


AlV 


7.0 


7.0 


DG Tau,RW Aur 




HD31464 


04 57 06.426 


+24 45 07.90 


45 


G5V 


8.6 


7.0 


DG Tau,RW Aur 




HD139364 


15 38 25.358 


-19 54 47.45 


53 


F3V 


6.7 


5.7 


HD 142666,HD 144432,AS 


205 A 


HD 141247 


15 48 11.770 


-04 47 09.684 


85 


F9V 


8.1 


6.7 


HD 141569 




HD 143459 


16 00 47.633 


-08 24 40.87 


141 


AOVs 


5.5 


5.5 


HD 142666,HD 144432,AS 


205 A 


HD 145788 


16 12 56.583 


-04 13 14.912 


171 


AlV 


6.3 


6.2 


HD 141569 




HD 149013 


16 32 38.133 


-15 59 15.12 


41 


F8V 


7.0 


5.7 


MWC 863 




HD 163955 


17 59 47.553 


-23 48 58.08 


134 


B9V 


4.7 


4.9 


MWC 275 




HD 170657 


18 31 18.960 


-18 54 31.72 


13 


KIV 


6.8 


4.7 


MWC 275 




HD 183324 


19 29 00.988 


+01 57 01.61 


59 


AOV 


5.8 


5.8 


V1295 Aql 




HD 188385 


19 54 40.200 


+07 08 25.27 


81 


A2V 


6.1 


6.1 


V1295 Aql 




HD199099 


20 53 26.390 


+42 24 36.72 


138 


AlV 


6.7 


6.7 


AS 442, MWC 1080 




HD1404 


00 18 19.657 


+36 47 06.81 


43 


A2V 


4.5 


4.5 


AS 442,MWC 1080 





References . — ( 1) [MuzeroUe et"aH ||2003|') : ('2') lwiiite fc Ghej toOlh : (S) lEisner et all to^*): (4) iMonnier et al] 1I2OO6I) : (5) 
lEisner et all l|2005i v Calibrator star distances are based on liirpa~pa;all ax measurem e^^ b^^^^liad ). 
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Table 2. 


Results of Modeling 














Source 




Tin 






Tin 






^dust 






^dust 






(K) 


(AU) 










(K) 






(AU) 


Dust Models 


RY Tau 


44.34 
















1750 ± 


10 


0, 


,16 


± 


0.01 


MWC 480 


31.85 
















1450 ± 


10 


0, 


,20 


± 


0.01 


MWC 758 


17.89 
















1610 ± 


10 


0, 


,17 


± 


0.01 


DG Tau 


4.70 
















1260 ± 


10 


0, 


18 


± 


0.01 


RW Aur 


3.60 
















1330 ± 


10 


0, 


,14 




0.01 


HD 142666 


1.55 
















1670 ± 


10 


0. 


,10 


± 


0.01 


HD 144432 


10.92 
















1420 ± 


10 


0. 


,20 


± 


0.01 


MWC 863A 


41.40 
















1360 ± 


10 


0, 


,28 


± 


0.01 


V1295 Aql 


25.20 
















1350 ± 


10 


0, 


,48 


± 


0.01 


AS 205A 


2.08 
















1850 ± 


10 


0, 


,14 


± 


0.01 


MWC 275 


19.93 
















1750 ± 


10 


0, 


,20 


± 


0.01 


AS 442 


8.63 
















1580 ± 


10 


0, 


,57 


it 


0.01 


MWC 1080 


24.39 
















2190 ± 


10 


1. 


,39 


± 


0.01 


Dust+Gas Continuum Models 


RY Tau 


0.30 


2684 ± 371 


0.02 ± 


0.04 


1.623 


± 


4, 


.655 


1105 ± 110 


0, 


,31 


± 


0.02 


MWC 480 


1.04 


3199 ± 391 


0.07 ± 


0.01 


0.081 


± 


0, 


,029 


1105 ± 9 


0, 


,28 


± 


0.01 


MWC 758 


1.73 


3462 ± 167 


0.07 ± 


0.01 


0.104 


± 


0, 


,013 


849 ± 53 


0, 


,35 


± 


0.03 


DG Tau 


1.48 


2188 ± 423 


0.01 ± 


0.01 


2.000 


± 


0, 


,001 


1137 ± 29 


0, 


,21 


± 


0.01 


RW Aur 


1.71 


7125 ± 981 


0.02 ± 


0.01 


0.009 


± 


0, 


,001 


1233 ± 26 


0, 


,15 


± 


0.01 


HD 142666 


0.13 


5577 ± 22197 


0.01 ± 


0.12 


0.175 


± 


0, 


,989 


1236 ± 736 


0, 


,14 


± 


0.07 


HD 144432 


0.54 


3910 ± 479 


0.07 ± 


0.01 


0.048 


± 


0, 


,015 


1031 ± 19 


0, 


,29 


± 


0.01 


MWC 863A 


0.67 


7714 ± 88 


0.07 ± 


0.01 


0.009 


± 


0, 


,001 


1128 ± 5 


0, 


,35 


lb 


0.01 


V1295 Aql 


0.58 


3115 ± 464 


0.14 ± 


0.01 


0.063 


± 


0, 


,029 


1106 ± 10 


0, 


,62 


± 


0.01 


AS 205A 


0.58 


6952 ± 3684 


0.01 ± 


0.01 


2.000 




0, 


,001 


1461 ± 68 


0, 


,18 


± 


0.01 


MWC 275 


0.13 


3983 ± 544 


0.06 ± 


0.01 


0.129 


± 


0. 


,045 


1230 ± 14 


0. 


,28 


± 


0.01 


AS 442 


2.70 


2521 ± 924 


0.29 ± 


0.04 


0.199 


± 


0, 


,282 


681 ± 1688 


1, 


,34 


± 


4.36 


MWC 1080 


0.67 


2790 ± 218 


0.50 ± 


0.01 


0.687 


± 


0, 


,173 


1271 ± 14 


2, 


,33 


± 


0.04 


Dust+Watcr Models 


RY Tau 


2.21 


5427 ± 654 


0.02 ± 


0.01 


0.238 


± 


0, 


,055 


1211 ± 4 


0, 


,26 


it 


0.09 


MWC 480 


2.46 


4659 ± 799 


0.02 ± 


0.01 


0.173 


± 


0, 


,084 


1246 ± 3 


0, 


,25 


it 


0.10 


MWC 758 


2.41 


4632 ± 738 


0.02 ± 


0.01 


0.201 


± 


0, 


,064 


1241 ± 4 


0, 


,24 


it 


0.08 


DG Tau 


1.56 


16907 ± 470 


0.01 ± 


0.06 


0.020 


± 


0, 


,273 


1204 ± 18 


0, 


,19 


it 


0.09 


RW Aur 


1.68 


6114 ± 2414 


0.03 ± 


0.01 


0.009 


± 


0, 


,000 


1226 ± 25 


0, 


15 


it 


0.01 


HD 142666 


0.25 


7334 ± 6398 


0.01 ± 


0.01 


0.757 


± 


1, 


,496 


1396 ± 5 


0, 


12 


± 


0.28 


HD 144432 


1.13 


6812 ± 3310 


0.01 ± 


0.01 


0.574 


± 


0, 


,711 


1225 ± 2 


0, 


,25 


± 


0.24 


MWC 863A 


2.55 


5090 ± 1149 


0.01 ± 


0.01 


0.187 




0, 


,103 


1218 ± 2 


0, 


,33 


it 


0.14 


V1295 Aql 


1.68 


5368 ± 1560 


0.02 ± 


0.01 


0.232 


± 


0. 


,053 


1222 ± 1 


0. 


,55 


it 


0.32 


AS 205A 


0.48 


4222 ± 2727 


0.01 ± 


0.01 


0.611 


± 


0, 


,602 


1631 ± 6 


0, 


,16 


it 


0.14 


MWC 275 


0.57 


5271 ± 897 


0.02 ± 


0.01 


0.246 




0, 


,059 


1409 ± 4 


0, 


,25 


it 


0.19 


AS 442 


2.92 


13212 ± 14239 


0.01 ± 


0.01 


2.093 


± 


7, 


,427 


1375 ± 3 


0, 


,67 


± 


2.54 


MWC 1080 


1.39 


10730 ± 2952 


0.06 ± 


0.04 


0.782 


± 


0, 


,750 


1600 ± 9 


1, 


,90 


it 


3.33 
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Table 3. Properties of Br7 Emission 



I Dust+Gas Model Dust+Water Model 

Source -fiBr7 -FBr7/-'^cont Rsri -FBr7/-Pcont 





(AU) 


(AU) 


RY Tau 


0.02 




0.04 


0.01 


± 


0.03 


0.07 


± 


0.01 


0.02 


± 


0.02 


MWC 480 


0.07 


± 


0.01 


0.08 


± 


0.02 


0.01 


± 


0.01 


0.09 


± 


0.03 


MWC 758 


0.02 


± 


0.01 


0.03 


± 


0.03 


0.01 


± 


0.16 


0.05 


± 


0.02 


DG Tau 


0.07 


± 


0.01 


0.16 


± 


0.04 


0.07 


± 


0.01 


0.16 


± 


1.13 


RW Aur 


0.07 




0.01 


0.06 


± 


0.04 


0.01 


± 


0.01 


0.06 


± 


0.05 


HD 142666 


0.05 


± 


0.01 


0.01 


± 


0.05 


0.05 


± 


0.01 


0.02 


± 


0.05 


HD 144432 


0.07 


± 


0.01 


0.04 


± 


0.03 


0.02 


± 


0.01 


0.05 


± 


0.05 


MWC 863A 


0.02 


± 


0.01 


0.04 


± 


0.01 


0.02 


± 


0.01 


0.05 


± 


0.02 


V1295 Aql 


0.01 


± 


0.12 


0.08 


± 


0.02 


0.06 


± 


0.02 


0.09 


± 


0.02 


AS 205A 


0.01 


± 


0.01 


0.01 


± 


0.01 


0.01 


± 


0.01 


0.02 


± 


0.04 


MWC 275 


0.01 




0.01 


0.05 


± 


0.02 


0.01 


± 


0.01 


0.06 


± 


0.02 


AS 442 


0.05 


± 


0.05 


0.04 


± 


0.19 


0.12 


± 


0.03 


0.03 


± 


0.06 


MWC 1080 


0.01 


± 


0.02 


0.05 


± 


0.02 


0.06 


± 


0.04 


0.04 


± 


0.02 



Note. — The uncertainties listed here are la statistical errors for our 
model fits. As discussed in i|4.1l these error bars are probably too small 
for i?Br7i although we can state confidently that -RBr7 < 0.1 AU for these 
sources. 



0.8 



0.8 



2.1 



2.2 

X,(|a,m) 



2.3 



2.1 



2.2 



2.3 



Fig. 1. — V"^ iWt) ^-iid fluxes {right) for several scans on the target V1295 Aql. Observations 
of a calibrator bracketing each target scan were used for calibration, and each line corresponds 
to one scan. The sohd curves show the results when HD 183324 was used as the calibrator, 
and dotted lines show results when HD 188385 was used. While we average all of these 
measurements together to determine the average V"^ and fluxes for a target, we plot the results 
for each scan and calibrator individually here to demonstrate statistical and systematic 
uncertainties in the data. Note that the strong feature at 2.165 /xm appears in all scans, and 
corresponds to Br7 emission intrinsic to the source. 
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1.00 
0.95 
0.90 
0.85 

0.80 
0.75 



HD1 67564 





HD171149(-0.1) 



2.1 
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Fig. 2. — Calibrated {^eft) and fluxes {right) for two unresolved check stars. The data 
for HD 171149 has been shifted down by 0.1 in both plots for ease of presentation. We also 
plot the fluxes predicted for these stars from Kurucz models (dotted histograms) . We applied 
the same calibrations to these stars as we did to V1295 Aql. However, these stars are > 30° 
away from their calibrator stars, substantially further than any of our targets are from their 
respective calibrator stars. The calibration uncertainties seen here may therefore be more 
severe than for any of our target data. 
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Fig. 3. — The fractional transmission of the atmosphere over Mauna Kea for 1.6 mm of 
precipitable water vapor and an airmass of 1.5. The curve has been smoothed to the spectral 
resolution of our observations. The strong absorptions on either side of 2.05 /xm are due to 
atmospheric CO2. 
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Fig. 4. — Calibrated V'^ (solid histograms) and fluxes (dotted histograms) for our sample. 
Channel-to-channel uncertainties are estimated to be 3% for V"^ and 5-10% for fluxes. For 
DG Tau and RW Aur, the noise in the measured V'^ is larger for A < 2.2 //m, owing to the 
low fluxes (< 1 Jy) at these wavelengths. 
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Fig. 5. — Fitted uniform disk diameter as a function of wavelength for our sample. The 
dotted line in each panel marks the wavelength of the Br7 transition of hydrogen. 
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Fig. 6. — Uniform disk diameter as a function of wavelength for RW Aur A {solid histogram) 
in the spectral region of the CO overtone bandheads. The measured fluxes in this region 
are also plotted {dashed histogram) . Dotted lines show the wavelengths of the CO overtone 
bandheads. The uncertainties on measured sizes and fluxes, which can be found in Figures 
ElandCl are ~ 5%. 
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Fig. 7. — The circumstellar components of measured V"^ and fluxes for our sample, as a 
function of wavelength, with the predictions of various models. Dashed curves show the 
predictions of models including a single-temperature ring of (dust) emission. Solid curves 
represent a model that includes a gaseous inner disk that emits continuum emission and Bt'J 
emission, in addition to the dust ring. Dotted curves show the predictions of models that 
include the effects of water vapor opacity in the gaseous emission. 
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Fig. [71 — continued. 
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Fig. [71 — continued. 
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Fig. 8. — Spectra of best-fit dust+water models, decomposed into dust {solid curves) and 
gas components [dotted curves). The parameters of the models are given in Tabled The 
gaseous component includes continuum emission from material hotter than 3000 K, water 
vapor emission from gas cooler than 3000 K, and Br7 emission. 



